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Abstract 

The  thermochemistry  of  simple  silicate  oligomers  in  aqueous 
sodium  silicate  solutions  is  rationalized  by  a  model  that  qualitatively 
predicts  equilibria  among  monomer,  dimer,  and  trimer  silicate 
structures.  Unlike  previous  models,  it  incorporates  the  influence  of 
both  silicate  bonding  and  ion-solvent  interactions  on  the  stability  of 
the  anions.  Although  semiempiricai  molecular  orbital  calculations 
indicate  that  silicate  bonds  become  less  stable  as  charge  increases, 
ion-solvent  interactions  tend  to  stabilize  highly  charged  and  compact 
silicate  oligomers.  Furthermore,  polymerization  is  favored  by  the 
formation  of  highly  stable  water  molecules  and  by  a  decrease  of 
charge  repulsion  within  the  silicate  structures.  Although  MNDO 
calculations  indicate  that  the  linear  structures  are  more  stable  than 
the  cyclic  structures,  ion-solvent  interactions  tend  to  stabilize  the 
cyclic  structures.  Polymerization  of  compact  silicate  oligomers  is  also 
favored  by  producing  more  water  molecules.  Implications  are 
summarized  for  larger  silicate  species. 
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1  Introduction 

The  crystal  structure  of  a  zeolite  is  highly  dependent  on  the 
composition  and  temperature  of  the  synthesis  solution.  Recent 
reviews  of  the  historical  development  and  mathematical  models  of 
zeolite  synthesis  emphasize  that  future  successes  in  the  design  of 
well-controlled,  continuous  synthesis  processes  will  depend  on  a 
detailed  and  quantitative  understanding  of  the  distribution  of 
solution  species.  Whether  zeolites  are  formed  directly  from 
dissolved  precursor  anions  [1-3]  or  from  rearrangement  of  the 
network  of  a  precursor  gel  [4,5],  the  population  of  the  various 
precursors  must  still  become  quantitatively  understood  if  zeolite 
synthesis  is  to  be  engineered. 

Using  29Si  NMR  experiments,  various  workers  have  obtain 
quantitative  distributions  of  silicon  among  more  than  20  silicate 
structures  [6-16].  It  has  been  shown  that  the  structure  of  these 
silicates  influences  the  nature  of  aluminosilicates  present  during 
zeolite  synthesis.  For  instance,  increasing  the  silicate  ratio  of  the 
initial  silicate  solution  increases  the  proportion  of  aluminum  present 
in  aluminosilicate  structures  during  zeolite  synthesis  [17]. 

However,  without  the  aid  of  a  physically  meaningful  model, 
NMR  measurements  alone  will  not  be  sufficient  to  fully  characterize 
the  solution  during  zeolite  synthesis.  29Si  NMR  spectra  are  relatively 
insensitive  to  the  degree  of  ionization  and  to  structure  beyond  the 
second  coordination  sphere.  Several  NMR  peaks  from  even  relatively 
simple  solutions  remain  unassigned.  Furthermore,  it  is  difficult  to 
collect  NMR  spectra  at  the  high  temperatures  characteristic  of 
hydrothermal  synthesis  [2,18].  Therefore  we  have  sought  a  rigorous 
thermodynamic  model  to  relate  the  concentration  of  various  silicate 
species  to  the  overall  solution  composition.  Previous  attempts  to 
develop  such  a  model  have  relied  on  empirical  assumptions  and  so 
are  unable  to  explain  why  specific  species  are  present.  Such  insight 
is  essential  to  construct  models  for  more  complex  species. 
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Although  silicate  and  aluminosilicate  solutions  of  technical 
interest  typically  contain  a  wide  range  of  oligomeric  species,  in  this 
work  we  begin  by  analyzing  simple  silicate  solutions  in  which  the 
only  possible  species  are  monomers  and  dimers.  In  fact,  even  such 
simple  solutions  can  have  up  to  12  different  anions!  Rationalizing 
the  behavior  of  even  monomers  and  dimers  can  yield  a  great  deal  of 
insight  into  the  reasons  for  construction  of  larger  oligomers. 

Although  several  workers  have  characterized  silicate  solutions 
[6-14]  and  have  even  discussed  possible  reaction  networks,  few  have 
addressed  how  to  explain  reaction  thermodynamics.  For  instance, 
Guth  and  Caullet  [19-23]  have  modeled  the  distribution  of  silicate 
anions  with  reasonable  success,  but  a  number  of  simplifications  were 
used.  For  instance,  the  polymerization  equilibrium  constants  were 
assumed  to  be  independent  of  the  nature  of  the  polymer  species.  At 
a  given  Si  site,  the  free  energy  of  ionization  was  assumed  to  increase 
linearly  for  each  successive  ionization.  Finally,  for  successive 
ionization  occurring  on  neighboring  Si  sites,  the  free  energy  of 
ionization  was  assumed  to  follow  an  empirical  algebraic  form  [23]. 
Although  this  model  accurately  accounts  for  the  distribution  of  small 
anions  in  dilute  alkaline  solution,  it  does  not  accurately  predict  the 
distribution  of  larger  species.  Moreover,  it  fails  to  explain  why 
specific  silicate  structures  (e.g.  certain  rings  and  cages)  are  formed. 
This  leads  us  to  reexamine  the  simplifying  assumptions  that  were 
made. 


One  weakness  of  the  model  may  lie  in  the  assumption  that  all 
siloxane  bonds  are  equally  stable.  De  Jong  et  al.  [30],  for  example, 
have  used  molecular  orbital  calculations  to  show  a  strong 
dependence  of  the  bond  stability  on  the  bond  angle.  Another 
weakness  lies  in  the  functionality  assumed  for  the  ionization 
equilibrium  constants.  As  will  be  discussed  below,  experimental  data 
that  the  free  energy  of  ionization  for  successive  ionization  does  not 
vary  in  a  linear  fashion. 
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In  another  effort,  Sjoberg  et  al. [25,26],  used  29Si  NMR  to 
characterize  silicate  species  in  fairly  complex  alkaline  silicate 
solutions  and  thereby  to  estimate  a  large  number  of  equilibrium 
constants  [27,28].  Although  these  constants  reasonably  reproduce 
some  experimental  results,  only  discrete  degrees  of  ionization  were 
assumed  to  be  present  for  each  silicate  species  in  the  solutions  even 
though  there  is  no  a  priori  reason  why  an  oligomer  should  not  exhibit 
all  of  the  possible  degrees  of  ionization.  It  would  be  valuable  to 
learn  what  species  are  actually  viable,  since  as  the  size  of  oligomers 
increases,  the  number  of  possible  species  (and  equilibrium  constants) 
increases  catastrophically. 

We  start  by  considering  the  simplest  alkaline  silicate  system, 
which  we  can  verify  by  NMR  to  have  only  monomers  and  dimers. 
Since  the  number  of  possible  species  is  small  (five  monomers  and 
seven  dimers),  the  number  of  independent  equilibrium  constants  is 
manageable  as  only  eleven  independent  chemical  reactions  are 
possible.  By  first  restricting  our  attention  to  these  few  reaction 
equilibrium  constants  rather  than  to  literally  thousands  that  would 
be  needed  for  complex  solutions,  we  both  ensure  accuracy  and  will 
be  able  to  observe  meaningful  trends. 

We  can  specify  the  eleven  independent  reactions  among 
monomers  and  dimers  according  to  the  following  scheme: 


M°  <->  M1'  <->  M2'  M}"  <->  M4' 

T 

i 

D°«-»D 1  ■«-*D2'<-»D3'<-»D4'  *->D5  <->D6 

where  Ml*  and  D‘*  represent  monomeric  and  dimeric  anions 
(balanced  equations  are  shown  in  reactions  1-3  below).  These 
reactions  are  depicted  systematically,  with  polymerization  and 
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ionization  being  considered  separately  so  that  meaningful  trends  can 
best  become  apparent.  We  do  not  suggest  that  this  reaction  scheme 
is  necessarily  kinetically  meaningful,  but  it  is  the  most  insightful  for 
thermodynamic  equilibrium  calculations. 

First  we  extract  AH,  AS  and  AG  data  for  monomers  and  dimers 
from  previous  work  which  we  carefully  select  so  as  to  avoid 
complications  resulting  from  reactions  among  too  many  anions.  We 
also  demonstrate  that  this  data  is  capable  of  predicting  the 
concentration  distribution  of  monomers  and  dimers  detected  directly 
with  29Si  NMR.  Most  importantly,  though,  we  rationalize  these 
thermochemical  trends  using  both  MNDO  and  ion-solvation  theory. 


2  Experiments  and  Calculations 
2.1  29Si  NMR 

Sodium  silicate  solutions  with  a  total  silica  concentration  of  0.2 
mol%  and  with  silicate  ratios  (  SiC>2  /  Na20  )  ranging  from  0.02  to  0.5 

were  prepared  from  high  purity  silica  gel  (Aldrich),  sodium 
hydroxide  (99.99%  pellets  from  Aldrich)  ,  and  distilled  and  deionized 
water.  29Si  NMR  experiments  were  performed  at  298°K  on  a  GE-500 
spectrometer  with  a  10  mm  sample  tribe.  A  paramagnetic  relaxation 
reagent  (Cr(HD)  was  added  to  all  solutions  to  reduce  the  29Si  spin 
lattice  relaxation  time  to  1-2.  The  Cr111  concentration  was  low 
enough  [8,9],  though,  that  there  was  no  significant  change  in  the 
anion  distribution.  All  chemical  shifts  use  Si(CH3)4  as  a  reference. 
Spectra  were  recorded  using  256  26  ps  (-90°)  pulses  with  a  pulse 

repetition  time  of  10s,  a  spectral  width  of  10  kHz,  and  a  line 
broadening  factor  of  5  Hz. 
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2.2  MNDO  calculations 

In  order  to  study  the  effects  of  silicate  bond  stability, 
calculations  were  performed  using  the  MNDO  Hamiltonian  introduced 
by  Dewar  and  others  [29-32].  MNDO  has  recently  been 
reparameterized  for  silicon,  so  it  is  useful  for  rationalizing 
thermochemical  trends  for  silicates  [33].  A  number  of  fairly 
successful  calculations  have  been  carried  out  for  molecules 
characterized  by  k  bonding  and/or  partial  ionic  bonding  such  as  Si02 
and  SiCl4  [33].  Though  the  expected  error  in  value  of  the  heat  of 
formation  is  about  33  kcal/mol  silicon  [32],  the  trends  suggested  by 
MNDO  calculations  of  silicate  anions  should  be  meaningful. 

It  is  important  to  bear  in  mind  that  the  stability  of  solvated 
silicate  anions  is  quite  different  from  that  of  anions  in  vacuum.  As 
an  example  of  how  important  this  proviso  is,  MNDO  has  been  used  to 
predict  a  stable  pentacoordinate  silicate  structure  [34]  even  though 
only  tetrahedral  structures  are  observed  experimentally  in  aqueous 
solution.  All  results  reported  in  this  paper  are  for  optimized 
structures  constrained  to  tetrahedral  silicon  coordination. 


3  Thermochemistry  of  monomers  and  dimers 

3 . 1  NMR  Measurements 

In  Figure  1,  the  spectra  obtained  for  dilute  solutions  with 
increasing  silicate  ratio  (decreasing  pH  value)  show  that  the  solutions 
contain  only  monomers  and  dimers.  The  distribution  was 
determined  by  integrating  the  spectra  and  is  shown  in  Figure  2  The 
dimer  is  first  formed  in  significant  concentration  at  a  silicate  ratio  of 
0.1,  and  thereafter  it  increases  in  concentration  as  the 
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silicate  ratio  increases.  The  test  of  whether  thermochemical  values 
from  the  literature  are  correct  will  lie  in  their  ability  to  reasonably 
predict  the  value  of  R  where  dimerization  begins. 


3.2  Selection  of  Thermochemical  Values 
3.2.1  Ionization  of  Monomers 

There  is  a  great  deal  of  discrepancy  among  reported  monomer 
ionization  constants  [25,26,35-40].  Some  workers  [36]  used  silicate 
solutions  which  are  too  concentrated  or  at  too  low  a  pH  to  ensure  that 
monomers  are  the  only  species  present  in  solution.  The  interference 
of  other  reactions  brings  the  thermochemical  values  into  serious 
question.  Other  workers  [25,26]  used  sufficiently  dilute  and  basic 
solutions,  but  their  calculation  of  ionization  constants  is  based  on  the 
assumption  that  only  a  few  of  the  five  degrees  of  ionization  are 
possible.  However,  Flint  et  al.  [37]  used  very  dilute  solutions  and  all 
monomer  anions  were  considered.  Moreover,  the  free  energies  of 
ionization  obtained  by  Flint  agree  with  those  reported  by  Mel’nik, 
Babushkin  and  Iler  et  al.  [39-41].  Therefore,  we  have  taken  Flint's 
values  as  the  most  reliable. 

The  standard  Gibbs  free  energy,  AG-  for  each  successive 
ionization  is  shown  in  Fig. 3a,  where  the  reaction  is  of  the  type 

Si(OH)mq°^  +  oh;  =  s,(OH)„.1q';;^  +  H,oaq  ( 1 ) 

While  successive  ionization  of  the  monomer  becomes  more  difficult, 
this  trend  levels  off  at  high  charge.  This  trend  runs  contrary  to  the 
assumption  made  in  previous  models  [19-23]  that  the  free  energy  of 
ionization  varies  linearly  with  the  degree  of  ionization. 

Enthalpies  of  reaction  have  also  been  found  from  variable 
temperature  investigations  [40],  and  Figure  3b  shows  that  the 
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enthalpy  of  successive  ionization  of  the  monomer  goes  through  a 
minimum.  The  entropic  contribution  to  the  ionization  of  monomers 
can  then  be  obtained.  Figure  3c  shows  that  the  first  ionization  is 
strongly  entropically  favored,  but  subsequent  ionizations  are  almost 
equally  opposed.  Notice  that  ionization  becomes  more  difficult  only 
because  of  the  unfavorable  entropy  trend,  but  the  free  energy  is 
modulated  by  the  enthalpy  trend. 


3.2.2  Formation  and  Ionization  of  Dimers 

Sjoberg  [27,28]  estimated  formation  free  energies  with  the  aid 
of  both  29Si  NMR  and  conductimetric  methods,  and  from  these  we  can 
extract  both  dimerization  and  dimer  ionization  equilibrium  constants 
corresponding  to  reactions  2  and  3. 

2Si(OH)4^  =  (HO)3SiOSi(OH)3a(7  +  H2Oaq  ( 2 ) 

Si2(0H)6.„0„:w  +  oh;  =  Si2(OH)5.nq‘"'a;  +  h2oo<  (3) 

Thus  we  find  from  Sjoberg's  data  that  the  free  energy  of  dimerization 
has  a  value  of  only  -0.7  kcal/per  mol  Si.  Figure  2-4  shows  the  free 
energy  of  successive  ionization  of  the  dimer.  The  dimer  ionization 
trend  is  similar  to  that  of  monomers,  but  the  dimer  is  more  acidic 
than  the  monomer  at  low  charge  and  then  becomes  less  able  to  ionize 
at  high  charge. 

Notice  that  that  free  energy  of  dimerization  (2)  is  negligible 
compared  to  the  free  energy  of  the  many  various  ionization  reactions 
of  types  (1)  and  (3).  This  indicates  that  the  distribution  of  monomers 
and  dimers  is  governed  primarily  by  the  stability  of  the  various  ions 
and  the  activity  of  the  hydroxide  ion.  Dimerization  occurs  only  when 
the  activity  of  the  hydroxide  anion  is  low  enough  that  the  hydroxide 
is  no  longer  driven  to  ionize  as  many  silanol  groups  as  can  only  be 
provided  by  monomers.  Therefore,  it  is  clear  that  we  need  to 
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understand  the  ionization  energy  in  great  detail  for  various 
oligomers. 


3.2.3  Confirmation  of  thermochemical  values 

In  figure  5,  we  show  the  calculated  composition  at  298°K  as  a 
function  of  silicate  ratio  at  a  fixed  Si02  concentration  (0.2mol%)  using 

the  data  from  figures  2-4.  It  is  assumed  in  these  calculations  that 
activity  coefficients  are  weak  functions  of  R  in  these  dilute  solutions. 
We  observe  reasonable  agreement  with  experiment,  but  the 
concentration  of  dimers  is  underestimated.  If  we  maintain  the 
ionization  constants,  but  increase  the  magnitude  of  AG  of  dimerization 
to  that  estimated  from  quartz  solubility  (-2.0  kcal/molSi  [19]), 
agreement  with  experiment  improves  (Fig.Sb). 

Another  possible  improvement  would  be  to  estimate  the  effect 
of  the  silicate  ratio  on  ion  activities,  but  so  long  as  the  ionic  strength 
of  the  solution  remains  low  enough  to  rely  on  the  extended  Debye- 
Hiickel  theory,  ion-ion  interaction  will  not  much  affect  the 
distribution  of  ions.  At  the  concentration  used,  we  expect  activity 
corrections  to  be  negligible  at  a  silicate  ratio  as  low  as  0.5  (where  the 
ionic  strength  of  the  solution  is  about  0.04M),  but  should  become 
increasingly  important  at  higher  R. 


4  Rationalization  of  Monomer  and  Dimer  Thermochemistry 

So  far,  we  have  only  selected  thermodynamic  values  for  simple 
species  from  experiments  and  shown  that  they  are  consistent  with 
NMR  for  dilute,  very  alkaline  solutions.  If  larger  species  were 
present,  many  more  thermochemical  values  would  have  to  be 
specified,  and  their  accuracy  and  precision  would  be  in  doubt 
without  some  understanding  of  the  behavior  of  even  the  simplest 
ammonia.  For  this  reason,  we  seek  to  rationalize  the  trends  in  AG,  AH, 
and  AS  using  MNDO  and  solvation  theory. 
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4.1  MNDO  Calculation  of  Ionization  Thermochemistry 

Figure  6  shows  the  heat  of  formation  of  linear  silicate  anions  in 
vacuum  as  calculated  by  MNDO.  The  values  are  reported  per  mol  Si 
rather  than  per  mol  anion  because  we  wish  to  compare  the  stability 
of  a  given  Si  unit  in  each  environment.  Each  anion  becomes  less 
stable  with  increasing  charge  largely  because  ionization  increases  the 
Coulombic  repulsion  energy  between  atomic  centers. 

Ideally  we  want  to  calculate  the  ionization  enthalpy  of 
monomers  in  aqueous  state,  (reaction  (1))  but  with  MNDO  we  can 
only  obtain  results  in  vacuum.  Estimated  enthalpies  of  ionization, 
calculated  from  data  in  Fig.  6,  are  shown  in  Figure  7  for  the  following 
type  of  ionization  reactions. 


s.(°H)mq,“^  +  oh;  =  si(OH)m.1q(:r^  +  H2oa?  (4  > 

where  subscript  vac  indicates  that  silicate  species  is  in  vacuum  state. 
The  reference  state  will  be  corrected  in  the  next  section. 

Though  MNDO  correctly  predicts  that  the  dimer  is  slightly  more 
acidic  for  the  first  ionization,  it  incorrectly  indicates  that  the  dimer  is 
more  acidic  at  high  charge  and  that  each  subsequent  ionization  would 
have  a  more  positive  enthalpy.  This  trend  conflicts  with 
experimental  observations  (cf.  Fig. 3).  This  discrepancy,  then, 
suggests  that  the  interaction  with  solvent  plays  an  important  role  in 
stability,  so  in  the  following  section  we  consider  how  to  estimate  the 
energy  of  ion-solvent  interactions. 


4.2  Ion-solvent  Interactions 
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To  study  solvation  effects,  the  extended  Born  model  of  ion- 
solvent  interaction  :s  used  [42]: 


AHi-h2o 


NAnZje0Hw  _  HAnzteQpw 
(rt  +  rw)2  (r,  +  rw)J 


NAUje0)2 
2{r^2rw)  * 


1  _ 

£w  Erf  dT 


where 


aHi-h2o  : 


"a 

Ms 

n 

Pw 


contribution  to  heat  of  formation  due  to  ion 

solvent  interaction 

radius  of  ion  and  solvent  molecules, 

respectively 

Avogadro  number 

charge  number  of  the  ion 

dipole  moment  of  solvent  molecules 

solvation  number  of  the  ion 
quadrupole  moment  of  water 


The  first  term  arises  from  the  interaction  between  the  ion  and  the 
dipoles  of  water  molecules,  the  second  term  from  the  ion  water- 
quadrupole  interaction,  and  the  last  term  (the  Born  contribution) 
from  the  work  of  transferring  an  ion  from  vacuum  into  a  continuous 
medium  with  a  dielectric  constant  ew.  It  should  be  noted  that  even 

the  extended  Born  model  is  at  best  semi-quantitative  when  applied 
to  complex  silicate  ions  in  aqueous  solutions.  Our  current  aim  is 
limited  to  qualitatively  rationalizing  the  thermochemical  properties 

of  the  simplest  silicates. 

As  a  first  approximation  we  estimate  the  silicate  anion  radius 
as  the  average  distance  from  the  center  of  a  silicate  structure  to 
oxygen  atoms  in  the  end  groups  of  the  structure.  Physically,  this 
means  that  we  assume  that  the  negative  charge  on  the  anion  can  be 

treated  as  if  it  resides  at  the  center  of  mass  of  the  anion.  The 

number  of  solvating  waters  is  initially  assumed  to  equal  the  charge 


of  the  anion.  In  appendix  A  we  discuss  more  realistic  assumptions. 
The  partial  derivative  in  the  third  term  is  negligibly  small  for  a 
dilute  aqueous  solution  at  298°K,  so  it  is  neglected. 

The  contribution  of  the  various  ion-solvent  interactions  is 
shown  in  Table  1.  Each  contribution  behaves  in  a  qualitatively 
similar  way,  though  the  Born  term  grows  more  quickly  than  other 
terms.  The  total  enthalpy  of  successive  ionization  of  monomers  and 
dimers  including  contributions  both  from  MNDO  calculations  and  ion- 
solvent  calculations  are  shown  in  Fig.  8.  The  trend  for  the  dimer  is 
now  qualitatively  similar  to  the  experimental  trend  in  Fig. 3b; 
subsequent  ionization  becomes  more  enthalpically  favorable  and 
then  reverses  at  high  charge.  This  is  the  consequence  of  competition 
between  unfavorable  anion  destabilization  (MNDO  prediction)  and 
favorable  ion-solvent  interactions.  However,  two  problems  remain: 
1)  this  calculation  predicts  that  ionization  of  the  monomer  becomes 
increasingly  favorable  without  reversing,  which  conflicts  with  the 
experimental  trend  (Fig.  3b),  and  2)  larger  silicates  are  more  difficult 
to  ionize  than  smaller  silicates,  which  is  also  incorrect. 

The  experimentally  derived  values  of  the  entropy  of  ionization 
(Fig.  3c)  reveals  clues  as  to  how  to  improve  the  calculation  of  the  ion- 
solvent  energy.  The  monomer  ionization  trend  shows  that  the  third 
ionization  introduces  no  more  order  to  surrounding  water  molecules 
than  does  the  second  ionization.  We  deduce  that  AS  is  governed  by 
two  opposing  trends: 

1)  Each  ionization  consumes  a  hydroxide  anion  and  produces  a  water 
molecule,  leading  to  a  favorable  entropy  change  of  19.3  cal/K.mol 

anion  (obtained  from  standard  values  of  H2O  and  OH  ).  This 

aq  aq 

entropy  contribution  is  constant  for  each  successive  ionization 
reaction. 

2)  Strong  ion-solvent  interactions,  on  the  other  hand,  cause  water 
molecules  to  be  ordered  in  the  solvation  sphere  of  the  anion.  Since 
each  ionization  produces  a  more  highly  charged  anion,  this  more 
charged  ion  sphere  interacts  more  strongly  with  surrounding  water 


molecules,  causing  the  water  molecules  to  become  more  ordered. 

For  this  reason,  each  ionization  tends  to  be  entropically  disfavored. 

Apparently,  though,  the  singly  charged  monomer  is  much  less  able  to 
order  water  molecules  than  the  hydroxide  anion  (Fig. 3c).  whereas 
multiply  charged  anions  order  water  better  then  does  the  hydroxide 
anion.  The  organization  of  water  molecules  around  a  charged  group 
should  not  be  appreciably  different  for  differently  sized  silicates,  so 
we  expect  that  the  entropic  behavior  of  ionization  of  dimer  to  follow 
the  same  trend  as  the  monomer. 

Fig.  3b  also  indicates  that  the  highly  charged  monomer 
interacts  no  more  strongly  with  water  molecules  than  does  the 
doubly  charged  monomer,  so  simply  increasing  the  solvation  number 
of  a  monomer  in  proportion  to  the  charge  number  provides,  of 
course,  a  poor  model.  In  fact,  only  a  few  water  molecules  closely 
surround  the  monomer,  whereas  other  water  molecules  must  interact 
from  a  greater  distance.  Because  the  ion-solvent  interaction  energy 
falls  off  rapidly  with  separation,  the  increase  of  the  monomer  charge 
will  not  proportionally  increase  the  degree  of  ion-water  interaction. 
Moreover,  the  extended  Born  model  treats  an  ion  as  a  structureless 
sphere;  this'  is  a  very  crude  assumption,  since  for  a  large  structure 
(e.g.  a  linear  trimer  or  tetramer),  it  is  more  reasonable  to  assume 
that  a  water  molecule  interacts  with  one  charged  silicon  tetrahedron 
in  the  structure  rather  than  with  the  whole  structure.  Details  are 
discussed  in  appendix  A,  where  we  illustrate  that  ion-solvent 
interaction  for  various  silicate  structures  is  strongly  model 
dependent.  The  current  model  fails  to  quantitatively  predict  the 
ionization  trends  mainly  because  it  has  few  parameters  to  model 
such  a  complex  interaction.  However,  it  is  clear  now  that  the 
stability  of  silciate  anions  arises  from  a  balance  of  bonding  versus 
ion-solvent  interactions. 


5  Rationalization  of  Linear  and  Cyclic  Trimer  Thermochemistry 


1  3 


As  the  solution  becomes  more  concentrated  and  less  alkaline, 
cyclic  silicate  species  are  formed  in  preference  to  linear  species.!7*1 U 
Figure  9  illustrates  MNDO  calculations  of  the  enthalpy  of 
polymerization  of  the  following  type: 


n  Si(OH)4vac  SinOn.1(OH)2n+2vac^Linear)  +  (n-l)H2Oaq  , 

and  shows  that  polymerization  is  generally  favorable.  According  to 
the  extended  Born  model,  ion-solvent  interactions  should  not 
influence  these  trends,  since  for  now  we  are  considering 
polymerization  of  only  uncharged  species,  but  more  subtle  ion- 
solvent  interactions  (such  as  dispersion  forces  etc.)  are  neglected 
here.  Polymerization  is  favored  mainly  because  it  produces  water 
molecules,  and  polymerization  of  larger  linear  oligomers  doesn't 
become  more  favorable  because  the  silicate  bond  becomes  less  stable 
for  a  larger  structures.  Polymerization  of  cyclic  structures  are  more 
favored  over  the  linear  structures,  because  they  produce  more  water 
molecules.  This  helps  to  explain  why  experiment  shows  that  linear 
polymerization  is  limited  in  solution  [6-12], 

While  uncharged  cyclic  structures  are  favored  over  uncharged 
linear  structures,  the  energy  of  ionization  is  much  larger  than  that  of 
polymerization,  so  it  is  not  yet  clear  why  cyclic  anions  are  favored 
over  linear  anions.  MNDO  calculations  of  the  enthalpy  of  ionization 
for  linear  and  cyclic  trimers  indicate  that  the  linear  trimer  ionizes 
more  easily  than  the  cyclic  trimer  [Fig.  10].  By  assuming  that  the 
entropic  contribution  to  ionization  is  similar  for  both  linear  and  cyclic 
trimer,  we  would  come  to  the  experimentally  incorrect  conclusion 
that  the  linear  silicate  oligomers  should  be  favored.  Therefore,  we 
must  again  consider  ion-solvent  interactions. 

Using  the  same  assumptions  as  used  for  Table  1,  one  can  show 
that  as  degree  of  ionization  increases,  the  stability  of  a  cyclic  trimer 
can  increase  more  rapidly  than  the  linear  trimer  (Fig.  11),  and  the 
cyclic  trimer  is  better  able  to  ionize  in  solutions  (Fig.  12).  Cyclic 
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silicate  oligomers  ionize  more  easily  because  the  ion  size  is  smaller 
than  for  linear  structures.  The  Born  energy  of  ion-solvent  interaction 
causes  it  to  become  more  stable  at  a  faster  rate  than  in  the  linear 
anion  as  charge  of  structure  accumulates. 


6  Conclusions 

Independently  derived  thermodynamical  values  accurately 
predict  speciation  of  monomers  and  dimers  but  do  not  follow  the 
simple  trends  that  have  been  assumed  in  previous  modeling  studies. 
We  found  that  MNDO  alone  could  not  predict  the  trend  in  ionization 
enthalpy,  but  inclusion  of  ion  solvent  interaction  explains  these 
trends  qualitatively.  We  are  still  limited  by  the  crudeness  of  the  ion 
solvent  model,  but  we  demonstrate  that  the  experimental  trends 
result  from  the  competition  between  bonding  and  ion  solvent 
contribution. 

The  trends  revealed  in  this  study  are  consistent  with  data 
reported  by  other  workers,  but  they  reveal  more  physical  insight 
regarding  silicate  speciation.  Four  general  rules  are  obtained: 

1 )  Each  subsequent  ionization  of  a  silicate  species  in  solution 
grows  more  difficult  in  a  complex  fashion  that  represents  the 
interplay  of  three  influences: 

a)  the  ordering  of  waters  of  hydration  (producing  a  regular 
and  predictable  AS  trend), 

b)  destabilization  of  the  molecules  as  charge  increases  (which 
can  be  predicted  with  MNDO  calculations),  and 

c)  favorable  energetic  interactions  between  the  ion  and  solvent 
(which  can  be  qualitatively  predicted  by  judicious  use  of 
extended  Bom  theory). 

2)  Ionization  is  enthalpically  more  favorable  for  larger  silicate 
species  when  they  are  only  slightly  charged,  but  it  becomes 
more  difficult  then  for  a  small  species  as  they  become  more 
highly  charged. 
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3)  Ionization  of  cyclic  species  is  more  favorable  than  the 
correspond  linear  species  due  to  shorter  range  stronger  ion- 
solvent  interactions. 

4)  Polymerization  is  favorable  because  it  produces  highly  stable 
water  molecules,  and  polymerization  of  cyclic  silicates  is  more 
favorable  than  the  linear  silicates  because  it  produces  more 
water  molecules.  The  thermochemistry  of  small  silicates, 
though,  seem  to  be  governed  by  both  bonding  stability  and  ion- 
solvent  interactions. 

Ion-solvent  interactions  are  very  difficult  to  model.  Future 
success  in  predicting  the  behavior  of  silicate  and  aluminosilicate 
species  will  rely  on  modeling  ion-solvent  interactions. 
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Appendix  A 


The  extended  Born  model  fails  to  quantitatively  predict 
ionization  trends  for  linear  silicate  oligomers  because  of  the  difficulty 
in  identifying  the  number  of  solvation  water  molecules  of  a  silicate 
anion  and  because  of  the  crude  treatment  of  a  silicate  anion  as  a 
charged  sphere.  To  address  these  problems,  we  can  modify  the 
calculations  in  the  following  way: 

1)  a  water  molecule  interacts  only  with  one  charged  silicate 
tetrahedron,  where  the  radius  of  the  tetrahedron  is  the  average  of 
the  four  Si-0  bonds, 

2)  the  number  of  water  molecules  that  interact  with  a  specific  silicate 
tetrahedron  increases  initially  in  proportion  to  the  charge  of  that 
tetrahedron,  but, 

3)  when  the  charge  grows  beyond  2,  the  solvation  number  increases 
more  slowly  (Fig.  13), 

4)  the  charge  distribution  on  the  dimer  is  optimized  to  provide 
minimum  AH^HjQ  (so  that  one  end  is  first  charged  up  to  3-,  then  the 

other  end  charged  to  it-). 

The  modified  calculations  (Fig.  14)  show  that  better  agreement 
with  experiment  can  be  achieved.  Fig.  14  also  indicates  that  large 
silicates  ionize  more  easily  than  the  smaller  silicates  when  they  carry 
low  charge,  but  large  silicates  are  harder  to  ionize  than  the  smaller 
silicate  when  highly  charged. 

However,  using  the  modified  ion-solvent  interaction  model 
derived  above,  the  calculated  values  for  the  enthalpy  of  ionization  do 
not  show  that  the  cyclic  trimer  ionizes  more  easily  than  the  linear 
trimer  [Fig.  15].  In  the  modified  model,  we  made  an  assumption  that 
a  water  molecule  interacts  with  only  one  specific  charged  Si  group 
and  neglected  the  interaction  with  other  Si  groups.  But  we  must 
count  on  this  interaction  when  we  want  to  compare  ionization  trends 
of  two  structurally  different  species.  This  is  illustrated  in  Figure  16. 
As  in  the  illustration  of  Fig.  16a,  a  water  molecule  can  have  strong 


17 


interactions  with  two  Si  groups  in  an  equally  charged  linear 
tetramer.  However,  in  Fig.  16b,  a  water  molecule  can  strongly 
interact  with  four  charged  Si  groups  simultaneously.  So,  ion-solvent 
interaction  of  cyclic  structure  is  stronger  than  for  the  correspond 
linear  structure. 

Inclusion  of  such  complex  interaction  in  the  ion-solvent 
interaction  model  overstretches  the  extended  Born  model. 
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Table  1  Sample  calculation  of  ion-solvent  interaction  enthalpy  for 
monomers  (kcal/mol  Si) _ 
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Figure  Captions 


1  2 9 Si  NMR  spectra  of  sodium  silicate  solutions  at  a 
constant  concentration  of  0.2mol%  Si02-  The  chemical 
shift  reference  is  the  monomer. 

2  Measured  distribution  of  monomer  and  dimer  in 
sodium  silicate  solutions  at  a  fixed  silicate 
concentration  of  0.2mol%. 

3  a  AG  of  successive  ionization  of  monomers  at  298°  K: 

(reactive  1).  Data  are  calculated  from  Flint  [37], 
Mel'nik  [39]  and  Babushkin  et  al  [40].  Error  bars 
reflect  scatter  of  value  reported  by  various  workers. 

3  b  AH  of  successive  ionization  of  monomers  (reaction  1). 

Data  calculated  from  Babushkin[40]. 

3c  AS  of  successive  ionization  of  monomers  (reaction  1) 
using  data  from  Babushkin[40]. 

4  AG  of  ionization  for  the  successive  ionization  of  dimers 
(reaction  3)  obtained  from  data  of  Sjoberg  et  al  [27,28]. 

5  Comparison  of  calculated  with  experimental 

compositions  for  silicate  solutions  of  constant  silicate 
concentration  of  0.2mol%.  (5a)  AGd=-0.7  kcal/mol  Si; 

(5b)  AGd=-2.0  kcal/mol  Si. 

6  Enthalpy  of  formation  from  elements  in  their  standard 
states  for  monomer  and  dimer,  obtained  from  MNDO 
calculation. 

7  Enthalpy  of  ionization  in  vacuum  (reaction  4) 
predicted  by  MNDO  calculation. 
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8  Successive  ionization  of  monomer,  dimer  in  solution, 
calculated  by  conbining  MNDO  and  extended  Born 
theory. 

9  Enthalpy  of  polymerization  from  uncharged  monomers 
to  uncharged  dimer,  linear  trimer  and  linear  tetramer 
in  vacuum. 

1 0  Enthalpy  of  ionizations  of  linear  trimer,  cyclic  trimer 
using  MNDO  calculation. 

1  1  Enthalpy  of  formation  for  isolated  (MNDO  calculation 
only)  and  solvated  (MNDO  +  ion-solvent  interaction) 
linear  and  cyclic  trimers. 

1  2  Enthalpy  of  ionization  for  the  linear  and  cyclic  trimer 
in  solutions. 

1 3  Modified  correlation  between  number  of  solvation 
water  around  a  specific  Si-tetrahedron  and  charge  on 
the  Si-tetrahedron. 

1 4  Successive  ionization  of  monomer,  dimer  in  solution, 
calculated  by  MNDO  and  a  modified  extended  Born 
theory. 

1  5  Enthalpy  of  ionization  for  the  successive  ionization  of 
linear  and  cyclic  trimer  in  solutions  by  MNDO  and  a 
modified  extended  Born  theory. 

16a  A  water  molecule  interacting  with  a  charged  linear 
tetramer  (Si407(0H^4) 
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16b  A  water  molecule  interacting  with  a  charged  cyclic 
tetramer  (Si4Oj(OH^  4  ) 
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Figure  3a 
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Figure  3b 
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Figure  16 
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